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Abstract. This paper studies analytically the stability of solitary waves in a generalized 
Boussinesq equation with quadratic-cubic nonlinearity. For general values of two parameters 
a and b determining the system, unstable waves may occur. If however, as in a situation for 
which this Boussinesq equation was recently proposed as a model for pulse propagation in 
nerves, (a, 6) belongs to a certain natural regime, then all possible waves are stable. 



1. Situation and results 

This note is directly prompted by the article [5] in which Heimburg and Jackson suggest the 
partial differential equation 

(1) v tt + (-v + av 2 + bv 3 ) xx + v xxxx = 

as a model for pulse propagation in biomembranes and nerves and argue that this model 
reflects certain properties of nerve axons better than the well known Hodgkin-Huxley and 
FitzHugh-Nagumo equations. Since the appearance of [5] this model has been studied in- 
tensely; for general aspects of these studies we refer the reader to the recent survey [Tj. The 
interest in equation (JTJ) rests on the fact that it admits solitary waves, i. e., traveling- wave 
solutions 

(2) v{x,t) = V{x - ct) with V(±oo) = 0. 

It is some of these solitary waves that Heimburg and Jackson propose as good representations 
for pulses in the abovementioned biological contexts. Now, as in order for this to be the case, 
the solitary waves should be dynamically stable, they and collaborators recently studied this 
issue computationally [10] and found that solitary waves are numerically stable in the case 
that the two parameters a and b occurring in (pQ) assume certain values that are significant 
for the concrete contexts they investigate. 

The present note gives a complete picture of the existence and stability of solitary waves in the 
extended Boussines^ equation ([1]) by analytical deduction. While the extreme cases a = 0, 
and 6 = have been well understood before (cf. [2]), no simple scaling argument applies to the 
case ab ^ 0. In fact, the literature does not seem to provide any concrete results concerning 
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the stability of solitary waves for generalized Boussines^ equations 
(3) v tt + (p(v)) xx + v xxxx = 

with p" non-monomial. Ours here rely on findings reported in [BJ. As [B], our argumentation 
follows Grillakis, Shatah, Strauss [3] and Bona and Sachs |2j in considering the so-called 
moment of instability's second derivative, the sign of which allows to conclude or preclude 
the existence of growing modes in the linearization of ([T]) around a solitary wave ([5]). 

We first characterize the set of all solitary waves that are possible for equation ([1]). 
Theorem 1. Equation ([!]) admits positive solitary waves of speed c if and only if a, b, c satisfy 



b > 0, a £ R and c 2 £ 
b < 0, a > and c 2 £ 



0,1), or 

max{0,l + ^},l) . 



It admits negative solitary wave of speed c if and only ifa,b,c, satisfy 

• b > 0, a £ R and c 2 £ 

• b < 0, a < and c 2 £ 

With 



0,1), or 

max{0,l + ^},l) . 




'(/:=-! h \ 1 — (1 — c 2 ) , and v :- 




each positive solitary wave has the respective value v as its maximum, and each negative 
solitary wave has v as its minimum. 

To give a precise definition of stability for this context, we write (pQ) as a system of first order 
in time: 

v± — u x = 0, 

( 4 ) ^ 

u t +p{v) x = -v xxx . 

Definition 1. [2\ A traveling wave (V, U) of (jl]) is called (orbitally) stable if for each e > 0, 
there exists a 5 > such that for any solution 

(v,u) £ (V,U) + C([0,T);H 3 (R) x H 2 (R)) 

of Q, closeness at initial time, 

\\(v,u)(;0)-(V,U)(-)\\ mxL 2<5 

implies, besides existence for all times (i. e., one may take T = oo), in particular orbital 
closeness at any time, 

mf\\(v,u)(;t)-(V,U)(- + o-)\\ H i xL 2 <e forallt>0. 
Definition 2. We call solitary waves of (JTJ Heimburg-Jackson pulses, if 

(5) b ^~\ a2 - 
The following is the main result of this paper. 
Theorem 2. All Heimburg-Jackson pulses are stable. 



2 We use this name in analogy with common terminology for the generalized Korteweg-de Vries equation 
v t + (p(v)) x + v xxx = (cf., e. g., [3]). 
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While there is no equivalence, for arbitrary generalized Boussinesq equations Q, between 
stability of constant states and stability of solitary waves (cf. [6], assertion (ii) of Theorem 
4a), the following seems enlightening for the family of equations under study. 

Theorem 3. Equation (P) is linearly wellposed at any constant state if and only if ([5]) holds. 

In other words, for ([T]), stability of constant states does imply stability of all solitary waves. 

We also show 

Theorem 4. (i) Assume that a > and 

(6) b > - 2 -a 2 . 

Then there are values < c* < c* < 1 such that while all positive waves of speeds with c 2 > c* 2 
are stable, all positive waves of speeds with c 2 < c 2 are unstable. Furthermore there are values 
< q, < cj < 1 such that all negative waves of speeds with c 2 < c? and all negative waves of 
speeds with c 2 > c 2 are unstable. 

(ii) Interchanging the roles of positive and negative waves, the same statement holds given 
© anda< 0. 

The transition, for positive waves, between stability for 'fast' waves and instability for 'slow' 
waves vaguely reminds of such a transition in the FitzHugh-Nagumo model, cf. [HUH]- 

Theorems 1 and 4 imply in particular that for certain choices of a and b violating ([5]), there 
are unstable solitary waves. 

Theorems 1, 2, 3, 4 will be demonstrated in Section 2. 

The following finding is useful for deciding (in-)stability of individual solitary waves for cases 
violating 

Theorem 5. For any solitary wave in (P) there is a simple algebraic expression 

fj,(a, b, c) 

depending only on the system parameters a and b and the wave's speed c such that the wave 
is stable [unstable] if fi(a,b,c) is positive [negative]. 

Section 3 comprises a proof of Theorem 5 and plots of \x that also illustrate Theorem 4. 

2. Proofs of Theorems 1 through 4 

As on the one hand the cases a = and 6 = are covered in the literature as mentioned 
above and on the other hand the transformation v — > —v is equivalent to replacing p(v) with 
—p(—v), we assume without loss of generality for the remainder of this paper that 

a > and 6/0. 

Proof of Theorem 1. With 

(7) F(,,c) = I(c 2 -l), 2 + V + ^, 
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Figure 1. Graph of F and level curves of I with p(v) = —v + av 2 + bv s for 
fixed speed c=|fora = |,6 = l (left), and a = 2, b = — | (right). 

a solitary wave satisfies the profile equation 

V" = -c 2 V - p(V) 
( 8 ) _ dFjy, c) m 

dv 

this equation admits the first integral 

I(V,V') = \v' 2 + F(V,c). 

In order for a solitary wave to be at least possible, (V, V) = (0, 0) must be a saddle point; 
this is the case if and only if 

c 2 < 1, 

which we henceforth assume. Theorem 1 follows directly (cf. Figure [1]) from the fact that 
besides at 0, F(., c) vanishes exactly at v and v. 

Solitary waves thus occur in families C V parametrized by their speed c. The key tool for 
stability considerations is the so-called moment of instability, 

/oo 
{ c V'f dx, 
-oo 

and our proofs of Theorems 2, 4, and 5 are based on the following fact. 
Lemma. The solitary wave C V is stable [unstable] if and only if the second derivative 

m"(c) 

of the moment at the respective speed c is positive [negative]. 

For this fact and the underlying theory, we refer the reader to O El HH [H] . 
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Proof of Theorem 2. Heimburg- Jackson pulses (with a > 0) are positive. As in [6], we 
obtain 

( c V') 2 dx = 2 / (-2F(v,c)) 1/2 dv 

-oo JO 

(-2F(v(c) - u> 2 , c)) 1/2 w dw, with iu := (w(c) - v) 1/2 . 

o 

Differentiating twice yields 

(9) m (c) = 2 / — — — '- dv 

Jo (-2F{v,c)) 6 / 2 

with 

F(v )C ) = \{c 2 -l)v 2 + a -^+ b -v\ 

F v (v,c) = (c 2 - l)v + aw 2 + bv 3 = c 2 v+p(v), 
F c (v,c) = cv 2 , 

U F v (v(c),c)- 

It is not difficult to verify that positivity of the integrand in ([9|) is equivalent to positivity of 

Q( v ) '■= T^ 3 + \ ay2 ~ v + c ^'( c ) - 1 + 7^ a v 
Now, one easily checks that 

Q(v(c))=0 and Q'(v) < - ^1 + for all v € R. 

This implies that Q is indeed positive on the interval (0,v(c)) and thus that m"{c) > 0. 

Proof of Theorem 3. Equation ([3]) is linearly wellposed at constant states vq G R if and 
only if every solution of the form w(x, t) = exp (At + iojx) , uj £ R, A G C, of its linearization 

wtt + p'(vo)w xx + w xxxx = 

has Re A < 0. Since for any such mode w, 

A 2 - p'(vq)uj 2 + w 4 = 0, 

this is characterized by p'(uo) = — 1 + 2afo + 36vq < 0. p' has no real zeros iff (0) holds. 

Proof of Theorem 4. Here, we have to consider two different cases. Consider first the case 
of a positive wave. Instability of standing waves and hence 'slow'ly traveling waves follows 
from the following observation: At c = 0, 

(10) m"(0) = 4 / - / '" ; „ /9 dv < 0, 



/ (-2F( V ,0))3/2 

since F(.,c) < in the intervall (0, v(c)). Continuity of integral and integrand implies then 
stability of waves with speed c 2 w 0; this observation is actually a special case of [7]. On 
the other hand, to prove stability of 'fast' waves, we apply Theorem 4 in [6]; translated into 
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the present situation, this theorem guarantees existence of a c* S (0, 1) such that all solitary 
waves of speed c 2 S (c* 2 , 1) are stable, provided that p'(0) < and p"(0) > 0; with 

p'(0) = -1 and p"(0) = 2a, 

this is obviously satisfied. 

Consider now the case of a negative solitary wave, i.e., b > and v(c) = min C V < 0. The 
considerations for c 2 < 1 slightly change as the moment of instability is now 

/oo /""(c) 
( c V'f dx = 2 (-2G(v,c)) 1/2 dv with now v(c) := -v(c) 
oo JO 



and 



G(v,c) = F(-v,c) = \{c 2 - l)v 2 - V + \v\ 



and its second derivative is ([9]) with F and its derivatives replaced by G and its derivatives. 
An obvious analogue of relation ([10]) keeps implying instability of waves with speed close to 
0. The following observation now shows instability of 'fast' waves with speed c 2 < 1. The 
quantities v = —v,m,m" extend to the limiting value c 2 = 1, with £>(±1) = (4a)/(36) and 
v'(±l) = -3/a. Thus, for all v £ (0,v(±l)), 

sgnm"(±l) = sgnv (2G(v, 1) (v + 2v(l)) - v (G v (v, l)v'(l) + G c (v, 1))) 

( -2a b 
= sgn ^— + -v 

< 0. 

Now by continuity of integral and integrand this implies m"(c) < for c 2 < 1. (Note in 
passing that in this case the minimum of the wave, and thus the wave's amplitude, do not 
tend to zero for c 2 — > 1). 

As can be seen in Figure [21 stable negative waves occur when k = b/a 2 is large enough. 

3. Proof of Theorem 5 and plots of sgn(m"(c)) 
We turn from estimating m"(c) to evaluating this quantity. Let 

g(k, c) := arcsin I - 



+ 2k(l 



g(k,c) := - In 

and 



+ 2k(l 



h(k, c) := 



-§ + V-2fc(l-c 2 ) 
4 (18/ec 2 -2 -9k) y/1 - c 2 



k (2 + 9fc(l - c 2 )) 
Proposition. If a > 0, the following holds with k = b/a 2 . 

(i) In the case ofb>0 and a positive wave, the assertion of Theorem 5 holds with 

4 [2 ( , 7T^ 



fi{a, b, c) = h{k, c) - — J - (g{k, c) - |) 
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(ii) In the case ofb > and a negative wave, the assertion of Theorem 5 holds with 

4 l~2 ( 7T\ 

ii{a, b, c) = h(k, c) - — W - [g(k, c) + -J . 
(Hi) In the case of b < and a positive wave, the assertion of Theorem 5 holds with 

4 r~2 

//(a, 6, c) = /i(/c, c) - —J--g(k, c). 

Remark. Note that with a > 0, there are no negative solitary waves in the case b < 0. We 
refrain from formulating the obvious analogue of Proposition 1 for the case a < 0. 



+ 




Figure 2. Plot of m" = in the case a, b > for positive (left) and negative 
(right) waves. The horizontal axis is k = b/o? and the vertical axis is c. The 
± signs refer to regions where m" ^ 0. 



Proof. As cases (ii) and (hi) can be treated analogously, we consider only case (i). By ele- 
mentary integration, we obtain 



m(c) = 2 
= 2 



v(c) 



(-2F{v,c)) 1/2 dv 



v(c) 







b 2 2 2 



1/2 



dv 



a 



26(1 - c 2 ) + ^a 2 

y 



+ 36 U j + 96 21 



6 

2x1/2 



arcsm 



-bv — |a 



>2 + 26(l- C 2)) 



After some slightly tedious calculations, this yields 



v=v(c) 



+ 



v=0 



m"(c) 



4 (I8fcc 2 - 2 - 9k) Vl -c 2 4 [2 



1 



fc (2 + 9fc(l - c 2 )) 



3fc V k 



— arcsm — 



3 U 



+ 2k(l - c 2 



7T 

2 



□ 
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